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CALIBRATION AND EVALUATION OF SKYLAB ALTIMETRY FOR 
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to 
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A. G. Mourad, Principal Investigator 
S. Gopalapillal , M. Kuhner and D. M. Fubara, Co-Investigators 

A. H. Byrns, Code TF6 - NASA/JSC Technical Monitor 


PROGRESS AND RESULTS 


During this period, the major activity consisted of preparation 
of a technical paper on "The Significance of the Skylab Altimeter Experiment 
Results and Potential Applications". The paper was presented by the Principal 
Investigator at the NASA Earth Resources Survey Symposium, held in Houston, 
Texas on June 8-12, 1975. A copy of this paper which will be published in 
the Symposium Proceedings is attached. 
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NEXT PERIOD 

During the next period, we plan to aubmit the Final report 
after receiving the approved draft from you. 

TRAVEL PLANS 


The Principal Investigator attended the NASA Earth Resources 
Sympos* j mentioned above. No plans for travel are expected during 
the ''.exw .^eriod. 
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TilE SICNIFICV.’CE OP THE SKYL/Ml ALTIMETER EXI’ERl?a:NT 
RESULTS AND POTENTIAL APPLICATIONS 

by A. G. Mournd, S. Cupnlapillni and M. Kuhncr 
Bnttcllc, Colunbus Laboratories , Columbus, Oliio 


ABSTRACT 


Tlic SUylab Altimeter Experiment has proven the capability of the altimeter for 
measurement of sea surface topography. The geometric determination of tlie gcold/mean sea 
level from satellite altimetry is a new approach having significant applications in many 
disciplines including geodesy and oceanography. A Generalized Least Squares Collocation 
Technique was developed for determination of the geoid from altimetry data. Tl>c technique 
solves for the altimetry geoid and determines one bias ten.i for the combined effect of sea 
state, orbit, tides geoid and instrument error using sparse ground truth data. The 
influence of errors in orbit and a priori geoid values are discussed. Although the Skylab 
altimeter instrument accuracy is about 11m, significant results have been obtained in 
identification of large geoidal features such as over the Puerto Rico trench. Comparison 
of tli 2 results of several passes sho\.»s that good agreement exists between the general 
slopes of the altimeter geoid and the ground truth, and that the altimeter appears to be 
capable of providing more details than are now available with best known geoids. The 
altimetry geoidal profiles show excellent correlations with bathymetry and gravity. 
Potential applications of altimetry results ts geodesy, oceanography, and geophysics are 
discussed . 


INTRODUCTION 

The Skylab S-193 altimeter experiment was the first of a series of altimeter experi- 
ments recommended by the "Will lamstown Study" \ and the NASA "Earth and Ocean Physics 
Applications Program" (EOPAP)s. The primary objective of the Skylab altimeter experiment 
was to determine the engineering feasibility of the altimeter instrument and demonstrate 
its capability for .neasurement of sea surface topography. 

Three manned Skylab missions--SL/2 , SL/3, and SL/4--provided data from the S-193 
system. Geodetic analysis of Skylab S-193 altimeter data from mission SL/2 EREP piss ?'4, 

6, 7, and 9 is the subject of this paper. Tlie overall objective of the Battelle investi- 
gation is to demonstrate the feasibility of and necessary conditions in using the altimeter 
data for the determination of the marine geoid (i.e., the geoid in ocean areas). Tlie geoid 
is the equipotential surface that would coinciae with "undisturbed" mean sea level of the 
earth's gravity field. "Undisturbed" is the condition that would exist if the oceans were 
acted on by the earth's force of gravity only and by no other forces such as those due to 

ocean currents, winds, tides, etc. Thus, determination of the gooid/mean sea level is 

basic to the understanding of the oceans and associated dynamic phenomena such as currents, 

tides, circulation patterns and, hence, air-sea interactions. Tiere exist many geoids 
which have been computed by various methods. Most of these lack the accuracy and quality 
required for many scientific and practical applications. 


1 


PrcBfnt methods of determining the p.eoid depend on the knowledge or mi jiMin ment of 
the d. toiled gravity field all over the earth and the use of satellite per turh.it i ons to 
describe the j;cneral field. Measurement of the p.ravity field all over the earth requires 
the use of land-based, shlphorne, airborne and ocean-bottom gravity instruments, a process 
uiiich will take decades to achieve an accuracy of the order of 1-3 meters. Present 
knowledge of the geoid on a world-wide basis is probably not better than 5-30 m. 

The use of altinctry for determining the geoid is, perhaps, most significant because 
of Its fundamental applications to geodesy, oceanography, geolt*gy, geophysics, navigation, 
national defense, environment, resource development and several other applications. 

Satelli'e altimetry offers the most expedient and accurate method for determining the gcuid 
independent of gravity measurements. An accuracy of ± 10 cm for determining the marine 
geoid from satellite altimetry Is the goal of the >iASA's E0PAP.‘ Tlic results obtained from 
the Skylab S-193 .itltretor experiment proved that the concept is feasible. Tnc extent of 
applications of titi.meiry data to various disciplines and use- will depend on the degree of 
accuracy achieved and correlation that can be made between these data and the parameters 
Involved in the application areas, 

CONCEPT OF GEOID DETEWllNATION FROM SATELLITE 
ALTIMETRY AND REQUIREMENTS 

Figure 1 shews schematic geocentric relations of the various surfaces associated with 
satellite altimetry, TM is the raw altimeter range which has to be corrected for laboratory 
instrumental calibration, electromagnetic effects, sea state, and periodic sea surface 
influences to give TS. S represents the non-periodic "sea level". CT and CE, the geo- 
centric radii of the altimeter and E, its subsatcllite point on the reference ellipsoid, 
are computed from satellite tracking information. EG is the bsolutc gcoidal undulation to 
be computed from this investigation, while SG is the quasi- stationary departure of the mean 
instantaneous sea surface from the geoid - the "undisturbed" mean sea level. It can be seen 
from Figure 1 that the required gcoidal undulations arc given by 

EG - ET - TM - M5 - SG (1) 

where, MS represents the sum of the calibration co.istants and the orbit uncertainties, if 
any. SG represents the deviation of the surface to which the measurement is made from the 
geoid. Since we do not have any information on SG which is not considered to vary sig- 
nificantly over the length of profiles corresponding to different submodes of observations, 
the sum (MS + SG) Is considered as the calibration constant. 

Hie basic requirements for determining the geoid from satellite sltimetry are: 

1. accurate orbit determination; 

2. precise altimeter instrumentation; 

3. ground truth verification data; 

A. methods of interpolating and extrapolating altimetry into 
unsurveyed are.'.s; and 

5, separation of the geoid from sea surface topographic effects. 

Hie geoid to be determined must be in absolute position or geocentric (i.e,, centered at 
the earth's center of mass) and have correct scale, shape and orientation in order to meet 
the goals of geodesy and also make contributions to the solution of problems in earth gravity 
modeling, geophysics, oceanography, etc. Correctness of stiapc depends on the precision of 
the altimeter and, in theory, absolute centering and orientation arc dependent on tlie 
satellite orbit ephcmcris. Tlic correctness of geoid scale requires that the orbit ephemeris 
and the altimeter either have no biases or systematic errors, or that such biases and 
systematic errors must be known to an accuracy better than the error tolerance of the geoid 
to be computed. Currently and for some time to come, those two scalar conditions cannot 
be met because of unknown systematic errors or biases in tracking station geocentric coordi- 
nates, the earth's gravity model, the tracking systems and tlie altimeter itself. TJiere is, 
therefore, a need for other sources of scale and orientation control, which can be satisfied 
by the use of good ground truth geoid data. Poor ground truth data will only result in a 
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coarse scale. Fine scale. If necessary, has to be introduced tl<rout;h some furm of marine 
geodetic control. Since there is no established marine geodetic control available at 
present, the geoid computed in this paper will only have a coarse scale provided by the a 
priori geoid height input. 

APPRQACl. AND MAT1!I:MATICAL MODEL 

The approach to this analysis consisted of two basic steps. TItc first step was to 
filter the noisy altimetry data using the Generalized Least Squares Collocation technique. 

In the second step, the filtered altimetry data were compared with the a priori ground truth 
geoid in order to determine the calibration constants and the altimetry geoid profile. The 
basic purpose of this comparison was to have the altimetry geoid data on the same scale as 
that of the a priori geoid. llic details of the procedure used in this analysis is xUscussed 
in Mourad, ct al.^ Tlie general mathematical model used in this analysis consists of the 
linear relationship between the residual altitude from the altimetry data, the a priori 
geoid, and the calibration constants. For the purpose of this investigation, only one 
constant term representing the cumulative effect of all the possible systematic errors 
associated with the residual altitude, is used as a calioration constant. 

Hie ba.^ic condition equation is 

D - R - N + AC + n^ (2) 

ill ^ 

where, R is the measured altimeter range, which is Intrinsically related to 

I 

1. the geocentric coordinates, X , Y , Z .of the satellite at the 

• I « 1 • I 

instant of measurement; 

2. the geoid undulation, N , referred to a given reference ellipsoid 

1 

at the subsatcllitc point; and 

3. the algebraic sum of the biases in all the measurement systems 
involved. 

Except for the noise term (n), equation (2) is the result of rearrangement of equation (1) 
with AC, which is the total bias term involved in the measurement, representi^'v sum 
(MS + SC) and with 

D - F (X , Y , Z , a, f) (3) 

1 lift (I 

representing the height of the satellite above the reference ellipsoid given as a function 
of the geocentric coordinates of the satellite at the Instant of observation and with a 
and f being the parameters (semi-major axis and flattening) defining the size and shape, 
respectively, of the assumed reference ellipsoid. It is important, however, that both 
D and N (Equation (2)] refer to the sai.:e reference ellipsoid. If they do not, they must 

1 I 

be made compatible by effecting the appropriate corrections as and when they arc necessary. 

The general theory of Least Squares Collocation is presented in Moritz,® and the 
details of its application to filtering the altimetry data arc described in Mourad, ct al^ . 
The filtered altimetry observations, s, arc given by: 

8 - N + AC - C (D - R) (4) 

• • 
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.wl)cre,C, • is the nuto-covariance matrix for the (^eoid umlulatirnH correspondint; to thu 
altimetry observations, 


: c - C,. + C,, (5) 

I 

1 with being the error covariance matrix for the altimetry observations. Equation (5) 
implies that the geoid undulations and the altimetry observations are stochastically 
independent. The subscript | is left out in equation (4) to indicate that the equation 
is in matrix notation and represents a set of several equations corresponding to equation 
( 2 ). 


The elements of the covariance matrix C,« arc evaluated from a table of numerical 
covariance function, from Tscherning and Rapp * , using linear interpolation mctliod. H>e 
covariance between any pair of undulations is given as a function of the spherical distance 
between these undulations. 

Tltc output resulting from filtering the altimetry data is a set of filtered residual 
altitudes given by 

8 - N + AC (6) 

If 8 and N are accurate, AC can be estiniated by evaluating equation (6) at any one point. 
However, in order to minimize the uncertainties, if any, in the ground truth geoid, this 
equation is evaluated at several points along the profile using Least Squares Adjustment 
principles, s is ccasldercd observable and AC and N arc parameters to be estimated. In 
matrix notation, the observation equation corresponding to equation (6) would be 

V + AV, + W - C (7) 

where, V is the vector of residuels on -s and V, is that on AC and N. A is the design 
matrix and W is the mis:losure vector obtained by evaluating (N + AC - s) with the 
'observed values' for s and the a priori ground truth data for N. The a priori values 
for AC can be assumed to be zero. If the weight matrices of the observables and the 
parameters arc P and Px respectively, the matrix solution for the residuals for the 
parameters is given by 

V, - -(A^PA + Pxj'^A^^PW (8) 

Then, the altimetry geoid undulations, N, are given by 

N - s - AC (9) 

If stochastic independence is assumed among observables and among parameters, P and P, 

will be diagonal matrices. Those elements of P corresponding to AC can be assumed to be zero. 

Preliminary examination of the Skylab altimetry data revealed that every time the 
submodc of the altimeter changed, there appeared to be a change in bias in the measured 
altitude. Consequently, separate bias terms were assumed for data from the segment of the 
profile observed in a single mode-submode combination, niese bias terms can bo recovered 
simultaneously by suitably modifying the design matrix.^ 

Consideration of separate bias terms for different segments may result in a discon- 
tinuity from one segment to the other. Hiis can bo rectified by constraining the closest 
ends of the two adjacent segments to have the same value of undulation even though these 
ends do not correspond to the same point. This is a reasonable assumption considering 
the accuracy of the Skylab altimetry system combined with the small variation in undula- 
tions over short distances of about 8-15 km. In this analysis, these constraints arc 
effected through sequential solution, wliich means that the effect of the constraints on 
the parameters is evaluated and added to the solution obtained from equation (8). Hie 
reason for this approach is tl>e convenience and efficiency resulting from the special 
structures oi the matrices involved. 




RF.SULTS AND CONCLUSIONS 


llie major results obtained can be divided broadly into two groups. One croup is 
concerned with the effects of errors inherent in the various input data such as tlie orbit 
ephemeris, a priori etc. Hie other consists of the results of the actual analysis 

of the data from the Skylab FREP passes Ob, 07, and 09, 

Results of Preliminary Data Analysis 

Hic results from the first group have been obtained from the analysis of some pre- 
liminary data from EREP pass 09 mode 5. The use of such preliminary data was necessary 
in the absence of actual data which wore not available at the time of the error analysis. 
Details of the preliminary analysis are described in Fubara and Mourad.^ Figures 2 and 3 
show only two important results. Figure 2 shows the influence of a priori geoid height 
on the geoid scale. It shows that such a priori geoid input and errors in them affect 
only the linear scale and not the shape of the geoid. This is due to the fact that any 
inherent error is modeled by a constant bias. The a priori geoid input was taken from 
Vincent and Marsh geoid’ (GC-73 in Figure 3). AA is the resultant altimetry gcoidal pro- 
file based on GG'*73 as a priori, input. Two different sets of errors were introduced into 
GG-73 to produce A-I and B-I. Thiti resulted in altrimetry profiles A-0 and B-0. It is 
obvious that AA (control experiment) is shape-wise identical tv A-0 and B-0. However, 
the scale of the calibration constant has changed considerably in the same altimetry 
profiles. Figure 3 shows that the scale and orientation of the computed altimetry geoid 
is highly dependent on the orbital data used. In this analysis, two sets of orbital data 
(computed independently using different methods) v;erc used. These resulted in the alti- 
metry gcoidal profiles J'J, and BE. The scale discrepancy was removed through the calibra- 
tion constant used. Howe\er, the difference in orientations remained. The close agree- 
ment between AA and GG-73 may be due to the use of basically the same gravity model and 
coordinate frame in their computations. 

To summarize, results demonstrated that: 

1. the precision of the altimeter ranges affect the shape of the geoid; 

2. the orbit uncertainties affect both the scale and orientation; and 

3. available a priori ground truth gcoids can give only a coarse scale 
to the altimet-^ geoid. 

Thus, to get an absolute geoid correct in scale and orientation, marine geodetic controls 
to offset t! p systematic errors in the orbit and altimeter data appear to be required. 


Results of Final Data Analysis 

The analysis of the altimetry data was made only for EREP passes 0^*, Ob, 07, and 09 
whose approximate locations in the North Atlantic Ocean are shown in Figure 4. The 
analysis was accomplished in three basic steps: 

1. filtering; 

2. estimation of the parameters; and 

3. graphical presentation of the results. 

The basic inputs for the first step are: 

1. the altimeter ranges and the exact correlated time of each measurement; 

2. the associated orbit ephemeris; 

3. the parameters of the reference ellipsoid; and 

4. the covariance function for the geoid undulations. 
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Hu* snlullitp altimetry data for tin* four pasiiei* were received on map.netlc lapeB 
from ItASA/JSC. TlieBe data con.nlst of el^ht a’limeter rauKv observations in framen at l.Oi 
second intcrvt.ls. However, for the purpose of this investit;at ion , the mean of the ei|,ht 
observations in each frame is considered as one observation. Tills assumption should not 
' egrade the results fur the following collective reasons; 

1. A frame of observations covers an effective area of about 6 km by 13 km, 
since the ground speed of the Skylab was about 7 km/sec and the radius 
of the radar foot print was about 3 km;° and 

2. Considering the accuracy of the altinetcr system on board the Skylab, 
the change in gcoid over an area of size 6 km by 13 km v*ould be 
insignificant . 

The Skylab Best Estimate Trajectory (SKYli. f) data are also available on a tape at inter- 
vals of exactly 1/8 of a second. Only the earth fixed geocentric coordinates of the 
Skylab and the time of obscrvatioi. arc input from these data. 

The best available estimates for the shape and size of the geodetic reference 
ellipsoid are given by: 

1. flattening - 1/298.255 

2. semi major axis diameter " 6,378,142.0 meters. 

Ihis is the same as the reference ellipsoid to which the Marsh-Vincent 1973 gcoid is 
referred.'’ Tlie covariance function for gcoid undulation is taken from Tscherning and Rapp.' 
Ihis is a numerical covariance function compatible with the reference ellipsoid chosen for 
the analysis. 

The residual altitudes which are filtered in the first step and the ground truth 
geoid undulations taken from the Marsh-Vincent geoio map^ form the input for the second 
step which is the estimation of the calibration constants and the geoidal parameters. In 
the third step, the estimated gcoid profiles arc plotted against time along with the 
ground truth profiles for easy comparison for .ihapc. 

The data in each mode of the altimeter wee observed in several submodes, each of 
which consists of several sub-submodcs of observations. Tne altitude measurement of the 
satellite rbovc the ocean which is the only g> odetic deta of interest, comes only from 
Modes 1, 3, and 5 of the Skylab altimeter Ihc me.. '.mum magnitude of the gcoid undula- 
tions, referred to the best available refere ice ellipsoid, is of the order of about 125-150 
meters. Therefore, any residual altitude c. more than a conscr\’ative estimate of 300 
meters, is an indication of instrument mal unction in the altimeter measuring system. 
Consequently, data were processed only fr those segments of the passes (f/4, , 07, and 09) 

corresponding to modes 1, 3, or 5 where .he absolute residual altitude is less than 300 
meters. It was found that the altimet' observations suitable for geodetic processing 
come from submodcs 0, 1, and 2 in mod s 1 and 5 and from submodes 3, 4, and 5 in mode 3. 

The results of the data analy is consist of: , 

1. a set of bias terms (Table 1) recovered for various segments of the 
four passes; and 

2. a set of geoidal profiles. 

Examples of these geoidal profiles and their corresponding bathymetry and gravity profiles 
are shovm in Figures 5 thrv 10. Hie geoidal profiles consist of: 
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1. ('round truth (dublicd linos) used in cstimntinR the bins 

terms and cstnblibhinn the scale; 

p,coid nroiiles computed from the unflltercd data usini; the 
omputpu bias terms (thin lines); and 

3. geoid profiles corresponding to the «.lltcred data. 

Since tlie altimetry data were fitted *o the ground truth geoid in the determination 
of the bias terms, any scale error in the ground truth Wuulu also result in error in the 
altimetry geoid. Consequently, any deviation of the altimetry geoid profiles from the 
ground truth profile will have short periods. In general, the deviation between the .wo 
sets of profiles is within about 2-3 meters with the following exceptions: In pass 
mode S (Figure 9) the deviation ranges from 0 to 12 m. In passes QU and 46 (Figures 6 
and /), the maximum difference is about 12 m in the Puerto Rico trench area. 

A close examination of these differences Indicatec that these extreme deviations 
occur in areas of special features such as trenches, ridges and sea mounts. Passes 
and 46 cross the Puerto Rico trench on the west side. Pass 47 mode 5 is along the 
western edge of the mid-Atlantic ridge while pass 47 mode 3 crosses the Puerto Rico trench 
at the Eastern end. 

These deviations iTuiy be due to several causes: 

1. Residual errors due to orbit uncertainties 

2. High fre ;ucncy component of the geoid not reflected in 
the ground truth data 

3. Possible (nadir) alignment errors which result in the 
departure of the sensor field of view from the nadir 

4. Influence of sea state, tide' and ocean circulation effects 

5. Possible inaccuracicc in th i ;omputation of the ground 
truth data 

6. Errors introduced as a result of scaling these data off 
small scale world maps. 

Most of the systematic errors caused by the above would be absorbed in the bias terms 
recovered from the data, especially because of the shortness of the segments for which 
separate bias terms were considered. The short periodic deviations caused by sea-state, 
tides, ocean circulation effects, etc., would be of the order of about 1-2 m. *'liis leads 
us to believe that, at least, the larger deviations are due to the short periodic compo- 
nents of the geoid not reflected in the ground truth geoid. This is confirmed by the 
fact that such deviations occur in the areas of significantly large geoldal features. 

Another striking difference noted in these profiles occurs in pass 47 (mode 5) which 
passes across the Puerto Rico trench (Figure 9). Hie trench, as indicated by both profiles, 
differs horizontally by al>out 30 seconds of time; equivalent to about 240 km. Looking at 
the gravity anomaly and bottom topography profiles, the altimetry geoid profile appears to 
be correct. However, further investigations may be needed in order to determine how a 
gross discrepancy such as this could have occured. 

Crossing of passes 44 and 46 (Figures 6 and 7) over the Puerto Rico trench and land 
mass area almost at the same place, has provided an ideal opportunity to compare the results 
to see the consistency of the altimeter system in determining the gold. The overlapping 
segments of the profiles for these passes are shown superimposed in Figure 8. Tlie agree- 
ment between the profiles is excellent. The only deviation at the beginning of pass 44 
mode 5, is due to the instrument transient response after switching pulse width, beam 
width and pointing between modes. Tliis aggrement indicates that the altimeter system is 
very stable and consistent. 
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From the rcRults nnd niinly*is prcnentcd thus far, somv p.^'Hornl con«m>nt s/ohsi-rval ions 
«nd conclusions can bo made: 

1. Hie procedure described and applied for filtering thu altimetry 
data and for estimating the calibration constants (bias) and 
geoid undulations has produced ve**/ satisfactory and realistic 
results; 

2. Ihc bias terms recovered for different segments of the same 

pass are significantly different. Tiicrc appears to be very little 
or no rorrelntion a‘< ''ng the bias terms associated with the same 
submodcs ; 

3. Hie agreement between the general slopes of the tvw geoid 
(altimetry and ground truth) profiles shows the viability of 
the altimetry tecliniquc in determining the marine geoid; 

4. Hie magnitude of the deviations of the altimetry geoid from 

the conventional ground truth leads to tlv- conclusion that these 
deviations are mostly due to the high frequency components 
of the geoid rather than due to other causes; 

5. Hie Skylab altimetry data analyzed here have provided ample 
evidence that the altimetry sensor is very sensitive to 
local geoidal features such as trenches, ridges, and sea 
mounts; 

6. Excellent agreement between the results obtained for the 
same place at different times shows that the satellite 
altimetry is precise and self consistent except for bias 
terms ; 

7. Hie correlation between the altimetry geoid and the gravity 
anomaly and the ocean bottom topography profiles have show. 

It to be useful in vcriiying major discrepancies in the 
conventional geoid. Hicse correlations may also be usefri 
in other applications such as in geology, geophysics, e'c. 

CORRELATION OF SKYUB ALTIMETRY WIHl 
BAHIYMETRY AND GRAVm 

Hiere are several obvious correlations of the Skylab altiroctrydctermincd geoid with 
ocean-bottom topography (bathymetry) and surface free air gravity anomalies. Each of 
Figures 5, 6, 7, 9, and 10 shows examples of these coi relations. In Figure 5, the 
altimeter responded clearly to the variation in topography and gravity anomalies caused 
by the Blake escarpment (about S m in geoidal heights, correlated with 3600 m depth changes 
and about 100 mgals in free air gravity anomalies). Note that the ground truth geoidal 
profile did not show such a change. The free air gravity anomalies show, also, some 
correlations with bottom topography over the continental shelf and the Blake plateau. 
Unfortunately, there are no gravity data (indicated by +-H- on the profiles) taken at the 
Blake escarpment. Figure 6 shows very strong correlations, over the Puerto Rico Trench, 
between the altimetry geoid and both the gravity and bathymetry, (about 13 m geoidal height 
change cor" nonding to about 6,000 m change in depth and 400 mgal in free air gravity 
anomalies). . ^ain, note the difference from the ground truth geoid. The altimeter also 
responded clearly to the land mass as evidenced by the sharp 'rise over Puerto Rico. 

Passes *4 and i‘6, over the Puerto Rico trench and Puerto Rico land mass, (Figures 6, 

7, and 8) provided an excellent opportunity to show the repeatability of the altimeter. The 

altimeter data also appear to show some changes due to the shallow water as the Skylab 
approaches the Lesser Antilles. There is also strong correlation between this altirctry 
geoid and the gravity, while the ground truth geoid is very smooth. Figure 9 (a segment 
of pass ^7) shows a strong correlation with the Puerto Rico trench (this is the eastern 
end of the trench wliich is different from those shown in passes (f4 and #6). Also the rise 
in the altimeter geoid correlates well with the rise in topography going toward the Lesser 
Antilles. Uov.’cver, the ground truth geoid appears to be shifted laterally by about 200- 
230 km. 
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imooth f,t'old over the whole prollle, xi»t* p,eoid variations do not conform to the I ot tom 
topocraphy ond pravity particularly over *he continental slope. It does correlate with 
the ground truth geold, however, with one exception where the ground truth profile sliowt 
m slight dip corresponding to the continental slope. Tlicrc is also a relative tilt between 
the altimetry gcoid and the ground truth gcoid. 

In sumnury, there are pronounced correlations bett^cen the altimeter-determined geoid, 
the topograpliy and gravity anomaly profiles, particularly those corresponding to special 
earth structures such as the Puerto Rico trench. In many cases, the altimeter appeared 
to sense many more details in the gcoldal surface than those reflected in the ground truth 
gcoid. Tnerc arc also some correlations with sea mounts or shallow near-surface topogro- 
phic features. McCoogan, et al showed similar correlations with the Puerto rico 

trench and the Mariana trench. In addition, he showed some correlations v;ith spcrlfl-. 
sea mounts and the Flemish cap as well as an overall correlation with the ground trw'.h 
geoid in around-thc-world-prof ile. Our investigation was limited to the data obtained 
from the four passes discussed in this paper. 

The causes of different gravity anomaly values, hence gcoidal features, associated 
with different but similar bottom topographic features are due to variations in the mass 
distribution in the earth crust and mantle and to the degree of isostatic compensation. 

In order to real' stically determine and evaluate the various deviations associated with 
the altimetry geoidal passes, more extensive effort is required. Su'*^ an effort could 
entail constructing sub-bottom profiles, based on various density assumptions and additional 
geophysical data (magnetic, seismic, geologic) and comparing the results with the gcoidal 
profiles. 


FUTURE ALTIMETRY APPLICATIONS 

As mentioned earlier, the extent of applications of satellite altimetry will depend, 
to a lorge extent, «n the degree of accuracy achieved in the determination of the geoid. 
With the CEOS-3 f ’.’'ite altimeter (launched April, 1975), a 1-5 ro geoid accuracy (abso- 
lute) Is expectr'. uy 1978/1979 another altimeter is being planned to fly on SEASA"^ with 
sub-meter expected accuracy in geoid determination. 

Future altimetry applications include the following: 

Ceodesv/n«v< cetion . - Any improvement in the determination of the geoid contributes 
to geedesy ond navigation. With an occuracy of ± 1 m in gcoid determination, many geodetic 
objectives can be achieved. The major contributions are in improving the determination 
of the size and shape of the earth and a unified datum and coordinate system on a world- 
wide '^asis. Determination of the figure of the earth, until the present, has depended 
largely on continental data. Satellite altimetry is certain to change that trend and 
include the ocean data covering over 70 percent of the earth's surface. The objectives 
of the now adjustment of the North American Datum (NAD) is 1 m. Tlie NAD is used for 
many civilian applications such as surveying, mapping, engineering operations, navigation 
and resource development. Tlie extension of control points and determining their three- 
dimensional coordinates to offshore areas as well as the determination of national and 
international marine boundaries must bo established in the same system. Accurate knowledge 
of absolute deflection of the vertical at sea, if combined with marine geodetic control, 
could provide the orientation required for all national datums, n^ese absolute deflections 
arc important for improving the accuracy of shipboard inertial navigation systems. Knowled 
of the geoid to ± 10 cm from altimetry could provide a worldwide reference lor the vertical 
datum which should contribute to investigations of land subsidence. 
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Cr.ivity field »lo t »M T'l n.u 1 on . - Hic r.rnviry field ',.:n hi* derived from n.ilcl li le nlli- 
mctry.JJ* A y.lobal t.oluti.Mi ui a ± I m nltinictiy i^cuid could resolve the y.ravity iield 
to about 1** X 1*^ whicli i ••preR.’nt s about 180 x 180 ^eopotential model. This resolution will 
contribute considerably lo inprovinc the p.eopotentlal model which is ubed in detennination 
of satcllitn orbit and tiissilc trajectories, llic altimetry yeoid cciild also provide infor- 
mation on the gravity aueinalier which in turn arc helpful, when combi-.ed with otiicr geo- 
physical data, in cxplor.it ion geophysics and identification of r,eoiogical structures. 

Mean s«'a level (M S..1. - a i 10 cm geoid will contribute signi f i.e'»nt ly to MSL determi- 
nation. Computation oi ilirection and magnitude of MSL slopes and the heiglits of each 
ocean relative to the otliors are important and still unresolved problems. Tliese problems 
arc further complicated because tht! results of geodetic and oceanographic computations of 
the parameters involved .ire different from each other. Determination of heights, directiorir 
and magnitude of sea level slopes relative to the continents is a key factor in land and 
environmental use and in siudies of the effect of changes in polar ice caps which affect 
marine life, meteoroli'gy and climate. 

plate tect o nics and ocean trenches . - Tliere is a correlation between the geoid/carth ' s 
gravity field a'’d gcophvslcal/geological phen'mena. Knowledge of these are required in 
earth and ocean physics studies. Some of these correlations wore shouTi in Skylab data 
above. To better understand and model this correlation, a knowledge of the fine structures 
of the ocean geoid is required. Studies of these phenomena contribute to better understanding 
of continental drift, polar motion and earthquake mechanisms. 

Occ inography . - With a ± 10 cm geoid. It is possible to determine the quasistationary 
departuiij of the sea surface topography from the geoid. is The departures of the sea 
surface topography from the geoid sudi as those due to tides, barometric pressure, wind, 
storm surges have practical applicatit i « in the determination of ocean dynamic phenomer.j 
such as circulation p.it terns, mass and nutrient transport, ocean tides, and ocean current 
influences. Most of tlicsc phenomena have important roles in monitoring and preserving the 
environment, in air-so.i interaction and in global numerical weather prediction. 

If the above parameters can be referenced to the geoid, satellite altimetry should 
contribute to the< r solution on a global basis. For example, sea slopes due to ocean 
current, could cau...e lee.il rise of water across the current on the order of 1 meter. 

Current slopes arc proportional to their speoj. Therefore, mass transport can be determined. 
The periodic effect of tides In the open oc jan is perhaps on the order of about 1 meter. 

Its determination is important particularly for the separation of the influence of earth 
tides. Also barometrie pressure could cause variations in the slope of the sea sunacc 
up to several meters. Other applications that cou’d be possible arc in the prediction 
of tsunamis (seismic sc.i wave) and storm surges. Tliese, however, may ,'e detected only 
If they occur during the pass of the altimetry Scatcllltes. Tsunamis amplitudes in the 
op'-n ocean range, pcrli.ips, from a few centimeters to about 1 meter. They have large 
wavelengths of the ordvr of several thousand km with about a 1-hour period. Storm surges, 
which is the local build-up of water due to distant violent storms (such as hurricanes 
and typhoons), could v.iusi' damage and reach wave heights of the order of several meters 
when they hit coastal .uvas. Their prediction and direction of movemenL could be of import- 
ance not only for coast.il .areas b t also for maritime ship operations. 

The various problems affecting the accuracy of achieving a ± 10 cm geoid must be 
solved in order to arrive at many of the above applications. The Skylab altimetry 
experiment, however, li.i.< demonstrated proof of concept of geoid determination from altimetry 
dala and that satellite .iltiractry is a potentially valuable tool having many useful 
applications. If Ci;0S-.« .ind SL^XSAT achieve their objectives, the Impact of altimetry on 
earth and sccan dynamie.s studies would be significant. 
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TABLE 1. “ RESULTS OF THE BIAS RECOVERY FROM ALTIMETRY DATA 


Pass Mode 

Suhmode 

Submode* 

Constrained 

Blns(tn) 

Filtered 

4 1 

0 

1 

-23.26 


1 

2 

-10.25 


2 

3 

-12.60 

5 

0 

7 

-28.12 


1 

8 

- 9.02 

6 5 

0 

1 

-24.03 


1 

2 

-11.63 


2 

3 

- 3.31 

7 1 

0 

1 

-52.43 


1 

2 

-40.28 


2 

3 

-79.28 

5 

0 

7 

-63.55 


1 

8 

-52.38 

< 

2 

9 

-46.56 

3 

3 

16 

- 


4 

17 

-08.87 


5 

18 

-87.39 

9 5 

0 

1 

-17.07 


1 

2 

- 6.07 


2 

3 

0.62 

3 

3 

10 

-19.99 


4 

11 

-20.53 


5 

12 

-18.95 


3 

18 

-36.95 


4 

19 

-36.26 


5 

20 

-34.80 


* Tlicsc numbers 
computat ionnl 

are asslpned 
convenience. 

sequentially in each 

for 
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Ficurc 2. - Effect of errors in a priori geoid height inputs and 
scale; dependency of calibration constant and gcoidal 
height on geodetic control (grount truth). 
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Figure 3. - Conventional Geoid and Satellite Altimetry Geoid Segments 

(Skylab SL-2 EREP Pass 9 Data) 
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Figure 5. - Gcoid undulations coniputed from Skylab pass-4 altinjctry 
data-mode 1. 
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Figure 6. - Geoid undulations -oaputed from Slylab pass-** 
altimetry data-modc 5. 
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Figure 7. - Geoid undulations coaputed from Skylab pass-6 
• altimetry data-mede 5. 
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Figure 9. 

- Geoid undulations computed 

from Sky lab pass- 7 
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Figure 10. - Geoid undulations computed from Skylab pass-9 
altimetry data-mode 5. 
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